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Resumen 


La investigación tuvo como objetivo ajustar una función de densidad de probabilidad (FDP) 
en plantaciones de Pinus caribaea var. caribaea Morelet Barret y Golfari (Pinaceae) de la Empresa 
Agroforestal Pinar del Río, Cuba. Los datos fueron tomados del Proyecto de Ordenación del Decenio 
2006-2016 de la Unidad Silvícola San Juan y Martínez. Para seleccionar los rodales incluidos en el 
estudio, se identificaron, según la Norma Ramal 595 y con la ayuda del SINFOMAP IV, aquellos 
con densidades superiores a 0,7 con diferentes edades y calidades de sitio. En total, fueron incluidos 
80 rodales pertenecientes a 41 lotes que tenían manejo recomendado. Se determinó la función de 
densidad de probabilidad de mejor ajuste con el software EasyFIT y como criterio para determinar 
el mejor ajuste se utilizó el estadístico de Anderson-Darling. La función de mejor desempeño resultó 
ser la de Weibull (2P) por mostrar mejor comportamiento en la predicción del número de individuos 
por clases diamétricas; las ecuaciones de los parámetros de escala y de forma se obtuvieron por 
regresión lineal múltiple por pasos. 


Palabras clave: función de densidad de probabilidad, raleos, redes neuronales artificiales, Pinus 
caribaea. 


Abstract 


The research had as objective to adjust a probability density function (PDF) in plantations of 
Pinus caribaea var. caribaea Morelet Barret and Golfari (Pinaceae) from Pinar del Río Agroforestry 
Company, Cuba. The data were taken from the Management Plan of the Decade 2006-2016 of 
the San Juan y Martínez Silviculture Unit. To select the stands included in the study, according 
to Branch Standard 595 and with the help of SINFOMAP IV, those with densities above 0.7 with 
different ages and site qualities were identified. In total, 80 stands from 41 plots with recommended 
management were included. The probability density function of the best fit was determined using 
the EasyFIT software and the Anderson-Darling statistic criterion was used to determine the best 
fit. The best performance function turned out to be that of Weibull (2P) because it showed a better 
behavior in the prediction of the number of individuáis by diametric classes; the equations of the 
scale and shape parameters were obtained by múltiple stepwise linear regression. 


Keywords: probability density function, thinning, artificial neural networks, Pinus caribaea. 


Introducción 

La caracterización de las distribuciones 
diamétricas proporciona información 
sobre la estructura del rodal y 
planeación de los tratamientos silvícolas en 
masas forestales bajo manejo. Conocer el 
número de árboles correspondiente a cada 
clase diamétrica de un rodal es de gran 
utilidad para el planeamiento y manejo 
de plantaciones, ya que permite estimar la 
variabilidad entre individuos, los posibles 
productos a obtener y, por lo tanto, el valor 
de la plantación. Las funciones de densidad 
de probabilidad (fdp) son utilizadas en el 
ámbito forestal para estimar la probabilidad 
de ocurrencia de un número de individuos 



y la probable distribución que presentará 
según sus diámetros, describiendo en 
términos matemáticos la estructura del rodal 
(Vallejos et al, 2004). Varios investigadores 
han realizado estudios en este tema entre 
ellos se encuentran: Nanang (1998) empleó 
la distribución Normal y Log-Normal, 
Zhang et al. (2001); Maldonado & Návar 
(2002); Gorgoso et al. (2007); Lei (2008), 

utilizaron la distribución Weibull, Álvarez 
& Ruiz (1997) emplearon la distribución 
gamma, y la distribución SB de Jhonson fue 
empleada por Kudus et al. (1999); Zhang et 
al. (2003) y Fidalgo et al. (2009). 

Pinus caribaea es una de las especies 
forestales que ha despertado un especial 
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interés a nivel mundial por su potencial de 
alta producción (Morales, 2011). En Cuba, 
es una especie priorizada en los planes 
de reforestación en la región occidental y 
central hasta el año 2020, por lo que se hace 
necesario una silvicultura efectiva que haga 
más eficaz su aprovechamiento (Castillo, 
2014). 

En las plantaciones de "pino" el 
incumplimiento de los planes de manejo 
provoca que los raleos no se ejecuten en 
el momento preciso, obteniéndose bajos 
rendimientos en los surtidos planificados 
según los objetivos de las plantaciones 
(Russo, 2015). 

La Unidad Silvícola San Juan y Martínez 
no está exenta de este deficiente control del 
raleo por lo que se hace necesario aplicar 
una herramienta que permita estimar el 
comportamiento de las clases diamétricas 
como herramienta para la toma de 
decisiones de las actividades silvícolas. 

La investigación se realizó con el objetivo 
de ajustar una función de probabilidad de 
densidad para la estimación del número 
de árboles por clase diamétricas en 
plantaciones de Pinus caríbaea var. caríbaea 
(Pinaceae) de la Unidad Silvícola San Juan 
y Martínez perteneciente a la Empresa 
Agroforestal Pinar del Río (Cuba). 

Material y métodos 

Se realizó un muestreo aleatorio 
estratificado, contando con el área basal 
como variable de estratificación. Se 
identificaron dos estratos: uno para las 
plantaciones que satisfacen la condición 
de un área basal y con una altura para 
densidad mayor que 0,7 y otro para los 
que no satisfacen esta condición, contando 
con una superficie total de 3 984,3 ha de las 
cuales 3111,3 ha son de áreas ocupadas por 
plantaciones con densidades menores que 


0,7 y 873 ha de plantaciones con densidades 
mayores que 0,7. Así, para el muestreo 
piloto, se tomaron 15 rodales con densidades 
menores a 0,7 y 10 para densidades mayores 
que 0,7 estableciéndose en cada rodal una 
parcela cuadrada de 400 m 2 , calculándose el 
número de muestra mediante la ecuación: 

n _ t 2 In=i WhSl 

E 2 + t 2 Y L WhSh 
b + r Lh=i N 

Fuente: Aldana, 2009 (1) 

Donde: 

n: número de muestras representativas 

Wh: proporción del estrato (h) en la 
población 

S 2 h : varianza por estrato 

N: número total potencial de unidades 
de la población 

E: error de muestreo. 

t: t de student 

Abordaje tradicional con las Funciones 
de Densidad de Probabilidades 

Para la determinación de la función de 
probabilidad de densidad de mejor ajuste se 
emplearon las redes neuronales artificiales 
recomendados por Bishop (1995) por tener 
un menor número de iteraciones debido 
a una mejor taza de convergencia y una 
capacidad de búsqueda más inteligente, se 
analizaron diferentes funciones (Tabla 1). 

Para el ajuste de los datos se empleó 
el software FitFD con el lenguaje de 
programación Jaba, uti izando como entorno 
de desarrollo IDE (Integrated Development 
Environment) Net Beans 7.1 y JDK 7.3 (Java 
Development Kit) disponibles de forma 
gratuita en el sitio Neuro Forest (http:// 
neuroforest.ucoz.com/); diseñado por 
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Tabla 1 . Función de Densidad de Probabilidades utilizadas para modelos de Distribución de Diámetros (MDD) 


FDP 

Modelo matemático 

Sb de 
Johnson 

/(*) - / 2níx y, g A2 l K+íi A«-J] [2] 

e, X, S, ey ^Parámetros del límite inferior (localización), amplitud 
(escala), asimetría (forma) y cutrosas 

Weibull (2P) 

f(X) = ^\-f 1 e~laf [3] 

a lai 

a = Parámetro de escala; /? = parámetro de forma; X = centro de 
clase de diámetro 

Burr (4P) 

roo- . y' *„ m 

4 i+ (tT] 

y y k parámetros de forma ( y > 0); /? = parámetro de escala (/? > 
0); a = parámetro de localización (a = 0) para distribución con 
dos parámetros 


Binoti et al, (2014). Para determinar el mejor 
ajuste se empleó la Prueba de Anderson- 
Darling. 

Resultados 

La caracterización de las distribuciones 
diamétricas proporciona información 
sobre la estructura del rodal y 
planeación de los tratamientos silvícolas en 
masas forestales bajo manejo (Quiñones et 

al, 2015). 

La función de densidad de probabilidad 
(FPD) Weibull es la más usada desde su 
introducción al campo forestal por Bailey 
& Dell (1973) y debido a su eficiencia y 
relativa facilidad de aplicación es una 
herramienta clásica de manejo en rodales 
puros coetáneos (Vanclay, 1994; Magaña 
et al, 2008). Su forma cerrada y gran 
flexibilidad es una ventaja porque puede 
adoptar diferentes formas desde una J 
invertida hasta distribuciones en forma de 
campana con diferentes grados de sesgo; 



su desventaja principal es su unimodalidad 
(Maldonado & Návar, 2002). Borders & 
Patterson (1990) generaron un método para 
caracterizar estructuras diamétricas de 
rodales que no requieren una distribución 
probabilística a priori. 

Los modelos de distribución de 
diámetros (MDD) o distribución diamétrica 
estiman el número de árboles por hectárea 
por clase diamétrica en edades presentes y 
futuras (Campos & Leite, 2013). Según Binoti 
et al (2014), la precisión de las estimaciones 
de clase de diámetro está influenciada por 
la elección de FDP que mejor describe los 
datos observados. 


Utilizando el sistema para ajustes de 
funciones de densidad de probabilidad 
software EasyFit se obtiene como la 
función de mejor ajuste la de Weibull de 
2 parámetros poseer el menor valor del 
estadístico de Anderson Darling (tabla 2). 
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Tabla 2. Bondad de ajuste de la Funciones analizadas y su ranking 


Función 

Anderson Darling 

Ranking 

Weibull (2p) 

0,32123 

1 

Johnson SB 

0,32429 

2 

Burr (4p) 

0,32458 

3 


Una representación de esta función para la muestra analizada se obtiene en la fig. 1 



Fig. 1. Ajuste de la función de Weibull de 2 parámetros para la especie en el área de estudio 


Con el empleo de la correlación las variables dasométricas obteniéndose la 
bivariada se correlaciona los parámetros matriz de correlación de Pearson (Tabla 3). 
de escala y de forma de la función de con 

Tabla 3. Matriz de correlación de Pearson 



| a 


D 

H 

Ho 

Do 

edad 

Dmin |Dmax 

a 

Correlación de 
Pearson 

i 

,567** 

-,402* 

-,356 

-,323 

-,435* 

-,103 

,938** 

,271 


Sig, (bilateral) 


,001 

,028 

,054 

,081 

,016 

,590 

,000 

,148 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

P 

Correlación de 
Pearson 

,567** 

1 

-,363* 

-,292 

-,212 

-,299 

-,270 

,796** 

,897** 


Sig, (bilateral) 1 

,001 


,049 

,118 

,260 

,109 

,149 

,000 

,000 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

D 

Correlación de 
Pearson 

1,402* 

-,363* 

1 

,915** 

,881** 

,802** 

,598** 

-,456* 

-,321 


Sig, (bilateral) 1 

,028 

,049 


,000 

,000 

,000 

,000 

,011 

,084 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

H 

Correlación de 
Pearson 

-,356 

-,292 

,915** 

1 

,965** 

,681** 

,690** 

-,398* 

-,252 


Sig, (bilateral) 1 

,054 

,118 

,000 


,000 

,000 

,000 

,029 

,178 


N | 

| 30 

30 

30 

30 

30 

30 

30 

30 

30 
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Ho 

Correlación de 
Pearson 

-,323 

-,212 

,881** 

,965** 

1 

,730** 

,693** 

-,335 

-,195 


Sig, (bilateral) 

,081 

,260 

,000 

,000 


,000 

,000 

,070 

,302 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

Do 

Correlación de 
Pearson 

-,435* 

-,299 

,802** 

,681** 

,730** 

1 

,499** 

,465** 

-,233 


Sig, (bilateral) 

,016 

,109 

,000 

,000 

,000 


,005 

,010 

,215 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

E 

Correlación de 
Pearson 

-,103 

-,270 

,598** 

,690** 

,693** 

,499** 

1 

-,221 

,470** 


Sig, (bilateral) 

,590 

,149 

,000 

,000 

,000 

,005 


,242 

,009 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

iDmin Correlación de 

1 Pearson 

,938** 

,796** 

* 

so 

lo 

r\ 

1 

-,398* 

-,335 

,465** 

-,221 

1 

,559** 


Sig, (bilateral) 

,000 

,000 

,011 

,029 

,070 

,010 

,242 


,001 


N 

30 

30 

30 

30 

30 

30 

30 

30 

30 

iDmax Correlación de 

1 Pearson 

,271 

,897** 

-,321 

-,252 

-,195 

-,233 

,470** 

,559** 

1 


Sig, (bilateral) 

,148 

,000 

,084 

,178 

,302 

,215 

,009 

,001 



N 

30 

30 

30 

30 

30 

30 

30 

30 

30 


**. La correlación es significativa en el nivel 0,01 (2 colas) 
*. La correlación es significativa en el nivel 0,05 (2 colas) 


En la tabla 4 se muestran los coeficientes 
de regresión parcial de las variables 
incluidas en el modelo de regresión, es decir 
la información necesaria para construir 
la ecuación de regresión en cada paso 
(incluyendo el término constante así como 
la significación estadística de los mismos) 

Tabla 4. Coeficientes del modelo 


rechazándose para el modelo 1 la constante 
siendo significativo todos los coeficientes 
para el modelo 2 con valores de tolerancia 
mayores que 0,01 lo cual rechaza la hipótesis 
nula para el supuesto de muticolinealidad 
como etapa de validación de modelo 
denotando la calidad del mismo. 


Modelo 

Coeficientes no 
estandarizados 

Coeficientes 

estandarizados 

t 

Sig. 

Estadísticas de 
colinealidad 

B 

Error 

estándar 

Beta 

T olerancia 

VIF 

1 

(Constante) 

,362 

,283 


1,281 

,211 



Dmin 

,324 

,023 

,938 

14,341 

,000 

1,000 

1,000 

2 

(Constante) 

2,653 

,272 


9,749 

,000 



Dmin 

,395 

,013 

1,144 

30,189 

,000 

,688 

1,454 

Dmax 

-,109 

,011 

-,368 

-9,713 

,000 

,688 

1,454 


a. Variable dependiente: parámetro escala ( a) 


Así resultan para la modelación del parámetro escala (a) como modelo de mejor ajuste el siguiente: 

a = 2,653 + 0,395 D min - 0,109 Dmax [18] 

En relación al ajuste del modelo del parámetro forma (fi ) se determinó por el mismo procedimiento 
que para el parámetro escala ( a) resultando la tabla 5. 
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Tabla 5. Bondad de ajuste del modelo para el parámetro forma (P) 


Modelo 

R 

R 2 

R 2 ajustado 

Error estándar de 
la estimación 

Durbin-Watson 

1 

,897 a 

,804 

,797 

2,10066 


2 

,965 b 

,930 

,925 

1,27460 


3 

o 

,953 

,948 

1,06353 


4 

,980 d 

,961 

,955 

,98823 

2,380 


a. Predictores: (Constante), Dmax 

b. Predictores: (Constante), Dmax, Drnin 

c. Predictores: (Constante), Dmax, Dmin, edad 

d. Predictores: (Constante), Dmax, Dmin, edad, H V 


Por su parte en la tabla 6 se muestran 
los coeficientes de regresión parcial de 
las 4 variables incluidas, para analizar la 
significancia estadística de los mismos, 
aceptándose para los 4 modelos la relación 


existente entre las variables al presentar 
sig.<0,05 lo que indica que el modelo 
mejora significativamente la predicción del 
parámetro forma. 


Tabla 6. Coeficientes no estandarizados para el modelo de parámetro forma (P). 


Modelo 

Coeficientes no 
estandarizados 

Coeficientes 

estandarizados 

t 

Sig. 

Estadísticas de 
colinealidad 

B 

Error 

estándar 

Beta 

Tolerancia 

VIF 

1 (Constante) 

,207 

1,967 


,105 

,917 



Dmax 

,725 

,068 

,897 

10,716 

,000 

1,000 

1,000 

2 (Constante) 

1,074 

1,200 


,895 

,379 



Dmax 

,532 

,050 

,657 

10,736 

,000 

,688 

1,454 

Dmin 

,403 

,058 

,429 

7,004 

,000 

,688 

1,454 

3 (Constante) 

-3,885 

1,711 


-2,271 

,032 



Dmax 

,602 

,046 

,744 

13,156 

,000 

,561 

1,781 

Dmin 

,394 

,048 

,418 

8,174 

,000 

,686 

1,458 

edad 

,106 

,030 

,172 

3,575 

,001 

,776 

1,288 

4 (Constante) 

-3,675 

1,592 


-2,308 

,030 



Dmax 

,640 

,046 

,792 

13,978 

,000 

,483 

2,070 

Dmin 

,338 

,051 

,359 

6,603 

,000 

,525 

1,905 

edad 

, 172 

,040 

,279 

4,285 

,000 

,365 

2,736 

H 

-,187 

,083 

-,142 

-2,261 

,033 

,396 

2,528 


a. Variable dependiente: B 


Así resultan sustituyendo los coeficientes de regresión no estandarizados el siguiente 

6 = - 3,675 + 0,64 Dmax + 0,338 Dmin + 0,172 edad - 0,187 H [19] 
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Discusión 


El análisis de las distribuciones 
diamétricas permite inferir el estado 
demográfico de las masas de bosque y sus 
posibles problemas de manejo. Los modelos 
de distribución de diámetros (MDD) o 
distribución diamétrica estiman el número 
de árboles por hectárea por clase diamétrica 
en edades presentes y futuras (Campos & 
Leite, 2013). 

Estos resultados del mejor ajuste de la 
función de Weibull se corresponden con los 
encontrados por Ferrere et al. (2001) para P. 
caribaea var. caribaea en Argentina así como 
para otras especies de Pinus como: P. pinaster 
por Álvarez & Ruiz (1998) en Galicia, España; 
P. pillea por Nanos (2002) en Valladolid, 
España; en plantaciones de P. durangensis, 
P. engelmannii, P. cooperi, y P. arizonica por 
Maldonado & Návar (2002) en Du rango 
México y para P. patilla por Santiago et al. 
(2014) en la región de Zacualtipán, Hidalgo, 
México; así como con los encontrados en 
plantaciones de P. taeda por Hirigoyen & 
Rachid (2014) en Uruguay, por Rubio et al. 
(2015) para parcelas permanentes de Pinus, 
Juniperus y Queráis en Nuevo León, México. 
Además, estos resultados discrepan con los 
encontrados para la especie por Montalvo et 
al. (1992) quienes encontraron como función 
de mejor ajuste la de Gauss. 


Como se puede constatar en la matriz 
de correlación existe una correlación fuerte 
y significativa del parámetro de escala 
con el diámetro mínimo y una correlación 
medianamente aceptable y negativa con el 
diámetro dominante y el diámetro medio de 
la plantación, por su parte para el parámetro 
de forma la variable más fuertemente 
correlacionada fue el diámetro máximo y 
también existe una correlación aceptable 
con el diámetro mínimo, estos resultados 
coinciden con los obtenidos para otras 


especies por Rennolls et al. (1985); Ortega 
(1989); Erviti, (1991); Maltamo et al. (1995); 
Alvarez (1997), Maldonado & Návar (2002). 

Con el modelo 4 se explica el 96,1% de 
la varianza del parámetro forma ( siendo el 
mejor. Al realizar el análisis o puntaje de la 
prueba de Durbin - Watson se observa que 
tiene un valor de 2,380 por lo que se acepta 
el supuesto de independencia serial de los 
residuos al encontrarse 1,5 y 2,5. 

El valor del nivel crítico permite afirmar 
que la incorporación de las variables edad y 
altura media contribuyen significativamente 
a explicar el comportamiento de la variable 
dependiente (parámetro forma ( forma. Las 
cuatro variables (Dmax, Dmin, edad, H) 
seleccionadas en el modelo final consiguen 
explicar el 96,1% de la variabilidad que 
presenta el parámetro forma (. 

Con una metodología similar a la 
utilizada para la obtención de los parámetros 
a y (3 por la relación de las variables de la 
masa en este trabajo, Rennolls et al. (1985) 
obtuvieron para Picea sitchensis, con datos 
de 120 parcelas, una varianza explicada del 
45,5% para un modelo lineal cuya variable 
independiente fue el dg, al igual Álvarez 
(1997) en P. pinaster en Galicia con R 2 
=0,99 y por García et al. (2002) en P. pinea 
en Valladolid, España quienes expresan la 
preponderancia de dg casi con exclusividad 
para explicar el parámetro a, así como para 
el parámetro (3 presenta un importante 
porcentaje de variación que no es absorbido 
por el modelo lineal que lo relaciona con dg 
y la edad. 

Kilkki et al. (1989) construyeron modelos 
de regresión lineal de los parámetros a y (3 
de la función de Weibull. En ambos casos, 
la variable de la masa que explicaba una 
mayor variación de los parámetros es dg en 
el caso del parámetro (3, los autores también 
incluyeron otras variables como g y la edad. 
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aunque no comentan la correlación entre el 
conjunto de las variables independientes, 
Ibáñez et al. (1995) obtuvieron un coeficiente 
de correlación de 0,79 en un modelo de 
predicción de a a partir de la raíz cuadrada 
del diámetro cuadrático medio 
logaritmo de la edad. 


y 


del 


Maldonado & Navar (2002) en 18 rodales 
plantados con P. durangensis, P. cooperi, P. 
engelmannii y P. arizonica obtuvieron los 
modelos del parámetro a en función de 
diámetro 1,30 medio y la altura dominante 
con un R 2 = 0,99, y el parámetro forma ((3) 
es explicado con un 92% del diámetro 
cuadrático basal, la altura media y el 
parámetro escala (a). 


de la EAF Pinar del Río, obteniéndose que 
para el parámetro de escala y de forma de 
la distribución por clase diamétricas fueron 
respectivamente: 

a = 2,653 + 0,395 Dmin - 0,\09Dmax 

/3 = -3,675 + 0,64Dmax + 0,338Dmin + 
0,172edad - 0,187 H 
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